clear all; close all; clc;
Rz=6400*10e+3; %радиус земли
H=19100*10e+3;%высота спутника над землей
R=Rz+H;%линия визирования
T=12*3600; %время полного оборота спутника вокруг земли в сек
dT=T/1000;%шаг по времени
t=0:dT:T;%массив временных отсчетов
fi=(t/T)*2*pi;%угловые отсчеты
x=R*cos(fi);%отсчеты координат
y=R*sin(fi);
Vx=diff(x)/dT;
Vy=diff(y)/dT;
V=sqrt(Vx.^2+Vy.^2);
Vtheor=2*pi*R/T;
e=V-Vtheor;
Rsv=[x;y];
Rr=[0;Rz];
for i=1:length(t)
    Rv(:, i)=Rsv(:, i)-Rr;
end

for k=1:length(t)
    Rvv(k)=sqrt(Rv(1,k)^2+Rv(2,k)^2);
end
Rvv(1001)=Rvv[];

%Vv=-(V*Rvv)/norm(Rvv);
